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Every completely positive map G that commutes which the Hamil- 
pH . tonian time evolution is an integral or sum over (densely defined) CP- 

maps Gfj where a is the energy that is transferred to or taken from 
the environment. If the spectrum is non-degenerated each is a 
dephasing channel followed by an energy shift. The dephasing is given 
• i-H , by the Hadamard product of the density operator with a (formally 

^ I defined) positive operator. The Kraus operator of the energy shift is 

■ a partial isometry which defines a translation on M with respect to a 

non-translation-invariant measure. 

As an example, I calculate this decomposition explicitly for the 
rotation invariant gaussian channel on a single mode. 

I address the question under what conditions a covariant channel 
destroys superpositions between mutually orthogonal states on the 
same orbit. For channels which allow mutually orthogonal output 
states on the same orbit, a lower bound on the quantum capacity is 
derived using the Fourier transform of the CP-map-valued measure 
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1 Introduction 



It was an important insight in the early days of quantum computing and 

quantum cryptography research that quantum theory should not be consid- 
ered as a purely physical theory, but it rather defines a new kind of informa- 
tion, called quantum information. The decisive feature is its fragility since 
quantum information is destroyed whenever one tries to copy it [1]. 

The question which physical processes and information channels preserve 
quantum information plays therefore a crucial role in the theory of quantum 
communication. A central concept for addressing this issue is the notion of 
a channel. Here we consider a channel or ("operation") to be a completely 
positive trace preserving map p i-^ G{p) where p and G{p) are positive op- 
erators of trace one acting on Hilbert spaces Hin and Tiout, respectively [2]. 
Analyzing quantum or classical channel capacities is in general a difficult 
task. Here we consider a certain kind of channels or operations, namely 
those which are time-covariant. For simplicity, we will assume that the input 
and the output state is a state of the same physical system. Then we have 
'Hin — 'Hout —■'H. Furthermore we have a Hamiltonian time evolution on Ti. 
generated by a densely defined self-adjoint operator H which reads 

a,(p):=e-^Ve^''*, (1) 

where p is a positive operator with trace one. We call a channel G time- 
covariant^ if it satisfies the equation 

atoG ^Goat. (2) 

It is a special case of the requirement 

G{U,pUl) = U,G{p)Ul , 

where g ^ Ug is an arbitrary unitary group representation. 

For compact groups, each covariant CP-map G has a representation (see 
[3] which refers partly to [4]) 

G{p)^Y.^,pL] 

3 

^Note that in the context of classical linear systems the corresponding property is 
usually referred to time- "invariance" . 
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where the Kraus operators Lj satisfy 

U,LjUl^J2dA9)Lk: (3) 

k 

and g i— > [djkig)]j,k<d is a d-dimensional unitary representation with arbitrary 
d (in [3] it is mainly focussed on the case that g ^ Ug \s an irreducible 
representation). Note that the time evolution can only be described by a 
compact group if it is periodic which cannot be guaranteed even in the finite 
dimensional case. In [5] the case is considered that G is replaced by a whole 
dynamical semi-group {Gt)te^+ which is covariant with respect to a group 
representation. 

Condition (2) appears naturally, for instance, in the following situations: 

1. Theory of timing information: Assume that the density matrix p 
is subjected to an arbitrary quantum operation G at a completely un- 
known time instant. Then the statistical description of this operation 
leads to a map G' which satisfies necessarily (2). This is the key idea in 
[6] where condition (2) defines a "quasi-order of clocks" which classifies 
systems with respect to their timing information. The time- covariant 
maps are exactly those which can be implemented without using ex- 
ternal clocks. Given a non-stationary state p, the set of states which 
can be obtained from p using time-covariant operations are those which 
have at most as much timing information as p. 

2. Decoherence: Dephasing of systems is described by a decay of the off- 
diagonal entries with respect to the energy eigenbasis. This channels 
satisfy clearly condition (2). 

3. Scattering processes / Quantum generalization of transfer func- 
tions: Time-covariant operations appear naturally in the description 

of scattering processes. Then G generalizes the scattering operator by 
including classical or quantum stochastic fluctuations. Here scatter- 
ing is understood in a rather general sense. Apart from the situation 
that a particle is scattered by the potential of another particle, one 
may, for instance, also think of a hght beam that passes a filter. One 
may consider time-covariant maps as quantum generalization of classi- 
cal time-invariant linear devices which are described by their transfer 
functions. 
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The paper is organized as follows. In Section 2 we recall briefly how to de- 
scribe classical time-invariant linear systems in signal processing by a transfer 
function [7]. To figure out to what extent this concept can be generalized to 
the quantum stochastic setting is the key motivation of this paper. In Section 

3 we recall abstract scattering theory in Hilbert spaces and explain in Section 

4 why timc-covariant channels are also considered as a natural generalization 
of scattering theory. In Section 5 we consider CP-maps which are given by 
the Hadamard product of the density matrix with a positive matrix. These 
maps will turn out as a building block of our decomposition. In Section 6 
we derive the general form of a time-covariant CP-map. The main result is 
that every time-covariant CP-map acting on a system with non- degenerate 
Hamiltonian has a (densely defined) decomposition as an integral or sum over 
(densely defined) CP-maps G^. is a dephasing followed by an energy shift. In 
Section 7 we address the question whether and to what extent covariant 
channels can destroy superpositions between states on the same orbit (with 
respect to the time-evolution). We use the general decomposition to derive 
a lower bound on the quantum capacity of a specific type of time-covariant 
channels, namely those which have as output mutually orthogonal density 
operators on the same orbit. In Section 8 we apply the decomposition to a 
single mode gaussian channel. 

2 Classical time-invariant linear systems 

To show the differences and the analogies to the theory of classical linear 
time-invariant channels we briefly rephrase the concept of transfer functions 
which is among the most important tools in classical signal processing. 

Consider a classical channel with an incoming and an outgoing signal. 
These may, for instance, be light pulses or electrical pulses, acoustical or 
other signals. Let t ^ Y{t) be the incoming signal where Y{t) is the value 
of an arbitrary physical quantity at the time instant t. Let Y'{t) denote the 
outgoing signal. In order to avoid problems with undefined Fourier trans- 
forms we assTimc that Y and Y' are vectors in L^(]R), the Hilbert space of 
square-intcgrablc functions on M. We assume that the system is described 
by a linear bounded operator A on L^(M) with 

Y'^AY. 
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Assuming timc-invariance of the channel we have Y'{t + s) = {AY){t + s) 
Defining the group of shifts {St)teR by 



{StY){s)^Y{t-s) 

we have 

StA ^ASt^teR. 

The imphcations of this condition can easily be derived by applying Fourier 
transformation to both sides of the equation: The shift operators St act then 
as the multiplication operators Mt 

Y^MtY 

with 

A linear operator commuting with all operators Mt is necessarily a multipli- 
cation operator as well (see [8], 12.1.5). Therefore A can be characterized by 
a so-called transfer function a with 

Y'{u) = a{u)Y{uj) . 

Note that it is essential that the quantity Y{t) is a scalar. If Y[t) is a vector 
of dimension d greater than 1 the spectrum of the time-evolution group St 
is degenerated and Y{uj) is an element of a vector space C"' or R"' and the 
transfer function a had to be replaced by a (i x d-matrix. It is straightforward 
to ask whether time-covariant quantum operations allow a natural general- 
ization of the transfer function. However, one should recall that the situation 
with non-degenerated spectrum is possible on the Hilbert space level (this 
is well-known in abstract scattering theory, as will be rephrased in the next 
section) but not on the level of density matrices, where the dynamics is gen- 
erated by the super-operator —i[H,.]. The latter has always degenerated 
spectrum. 

3 Abstract scattering theory 

A similar approach as above applies to a scattering process of a quantum 
mechanical particle: A particle which comes from the infinity and passes a 
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scattering potential. After it has left the potential it disappears to infinity. 
Here we only recall some standard results of scattering theory [9, 8] which 
are essential for this paper. Let H — + Hq he the total Hamiltonian of 
the system which consists of the free Hamiltonian Hq and the interactions 
term Hi. One assumes that the particle moves approximatively according to 
its free Hamiltonian Hq for t ±00. For potentials which decay sufficiently 
with the distance between particle and the scattering center, the limits 

lijj^) := lim exp(iHt) exp(—iHot)\'ilj) 

t— >±oo 

exist [9] on an appropriate subspace of wave functions There is a rich 
literature addressing the question under which circumstances and on which 
subspaces one can define a unitary scattering operator^ 

However, if it exists it commutes with the free Hamiltonian evolution 
exp(— ii^o^)- Therefore the free time evolution and the scattering operator S 
can simultaneously be described by multiplication operators. The analogy to 
Section 2 can even be closer in a scattering process with so-called Lax-Phillips 
evolution [8]. Assume that there exist subspaces Ti- < 7i and H+ < 7i 
( "incoming and outgoing subspaces" , respectively) such that 

for all negative t (the particle comes from the infinity) and 

e-'"^n+ c n+ 

for all positive t (the particle disappears to the infinity) . Assume furthermore 
that the intersection 

vanishes and that the span of all spaces cxp(— ii7t)?i+ as well as the span 
of all exp{—iHt)T-C- is dense in Ti.. Then one can assume without loss of 
generality that 7i is the set of square integrable functions on M with respect 
to the Lebesgue measure, and H is the multiplication operator {Hip){uj) — 
uj'4){uj). The scattering operator is then (like the transfer function in Section 
2) given by {Sip){u!) — s{u!)ip{u!) with an appropriate function s. 

^It is usual to denote the scattering operator with S. On the other hand, it is usual to 
denote shifts by St with some index t. He hope that this will not lead to confusions. 
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4 Quantum channels from scattering processes 



The unitary scattering operator S in Section 3 defines a channel G{p) — SpS''' 
which satisfies obviously the time-covariance condition (2) with respect to the 
free evolution at{p) = exp{—iHot) p exp{iHot) due to SHq = HoS. It is clear 
that scattering with an unknown potential could lead to an operation which 
is a statistical mixture of maps SpS^. But not only classical fiuctuations of 
the scattering potential lead to a mixture of output wave functions; due to 
quantum fiuctuations of the potential one may also have a CP-map which is 
not a mixture of unitary scattering processes. To see this, consider a bipartite 
quantum system with Hilbert space Ha^'Hb that evolves according to a joint 
Hamiltonian Hab- Assume that the joint evolution coincides asymptotically 
with the separate evolution generated by Ha + Hb and that the limits 

:= hm cxp{itHAB) exp{—iHAt) exp{—iHBt)\ip) 

exist in an appropriate sense. Then we have formally the same situation as 
in Section 3 with the substitution 

Hq Ha + Hb , H i-^ Hab ■ 

If a scattering operator U on Ha 'SiHb exists we have 

[U,Ha + Hb]^0 (4) 

Let Pa be an arbitrary initial density matrix of system A. Then the scattering 
process defines a completely positive map on the system B by 

G{pb) :=trA{U{pA0PB)U^). 

Due to (4) it is easy to verify that G is time-covariant if the initial state pA 
is stationary with respect to the time evolution exp(iif^i). Here we do not 
address the difficult question in which situation the hmits above exist. The 
problem of the existence of scattering operators in the quantum stochastic 
setting is addressed in [10]. This will not be our subject. The remarks above 
were only to show that time-covariant channels appear naturally also in the 
description of (possibly inelastic) scattering processes. Another situation 
where time-covariance appears is when the energy spectrum is discrete and a 
weak interaction with the environment is switched on. Then the interaction 
implements clearly a unitary U on system -|- environment such that U satisfies 
(4). 
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5 Quantum channels from Hadamard-products 



Here we consider a simple type of time-covariant CP-maps which will play a 
crucial role in the description of general time-covariant CP-maps. 

The Hadamard product A*B of two nxn- matrices A, B is defined as the 
entry-wise product {A* B)ij := AijBij. Remarkably, the Hadamard product 
between a density matrix p with any positive matrix M of the same size n 
defines a completely positive map: 

If M := ^ spectral decomposition of M then we may 

define diagonal matrices Dj which have the coefficients of the vector dj € 
as diagonal entries. Then the map 

p^DjpD] = {\dj){dj\)*p 

is obviously completely positive. So is the positive linear combination of 
those maps. Channels of this type have already be considered in [11] to de- 
scribe decoherence. It is clear that they commute with the Hamiltonian time 
evolution if the Hadamard product is calculated with respect to the energy 
basis. In the finite dimensional case, it we obtain also time-covariant chan- 
nels by the following construction: Let E := {cui, . . . , cUn} be the eigenvalues 
of H acting on C". Let \u}j) be the corresponding eigenstates. For any cr G M 
let ji, . . . ,jk be all indices j such that cuj -|- cr is in the spectrum E. Define 
the "partial shift" by 

■='^\uJj^+(T){uJj^\, 

where the sum runs over all spectral values uj with c<j + o" G E . 

The map p i— > S^rpS^ is time-covariant and also every map of the form 

where the sum runs over all a for which a non-zero partial shift exists. Each 
Mg. is an arbitrary positive matrix. It is easy to check that this map is 
trace-preserving if and only if 

^M^{u},iu) = 1 Vo; e E, 

a 

where the sum runs over all values a for which there exist an a; G E such 
that a; -|- cr G E. 
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In the following we will show that this is the most general form of time- 
covariant CP-maps. As will be shown, this holds in principle even for the 
infinite dimensional case with continuous spectrum when the sum is replaced 
by an integral over a potentially uncountable number of energy shifts. 

6 The general form of time-covariant chan- 
nels 

To understand our construction it is helpful to consider the finite dimensional 
case first. Using the Kraus representation [2] 

G{p) = J2 ^^P^] ' (5) 

3 

one may choose Aj such that they are eigenvectors of at (Since the derivation 
for the finite dimensional case follows actually from the general derivation 
later in this section we only mention briefiy that this can be derived as fol- 
lows. One shows that the representations t ^ [dkj{ty\ of (M, +) in cq. (3) 
corresponding to time-translations can be choosen such that they are irre- 
ducible, i.e., one-dimensional). The eigenvalue is exp(— iajt) where Uj may 
be any possible frequency difference uj — uj' . For each possible value cr we 
obtain a CP-map 

Let Aj = Uj\Aj \ be the polar decomposition of Aj where the partial isometry 
Uj is computed from the pseudo- inverse by 

One has O-tiUj) = exp{—iajt)Uj and equivalently 

[H, Uj] = ajUj and [H, \Aj\] ^ (6) 

since \ Aj \ is at-invariant according to at{Aj)at{Aj) = e~^'^^^e^'^^^. Using equa- 
tion (6) one checks easily that 

Ujluj) = Vj\u! + a) Va; e S 
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where Vj is an appropriate diagonal operator. Hence we can write 

with the diagonal matrix Dj :— \Aj\Vj. Using the remarks of Section 5 we 
may write 

G^ip)=S^iM^*p)Sl 

where = \dj){dj\ is defined as in Section 5 from the vectors \dj) of 
diagonal entries of Dj. 
We conclude: 

Theorem 1 (Decomposition in finite dimensions) 

Let G be a CP-map on which commutes with the time evolution at corre- 
sponding to a non- degenerate diagonal Hamiltonian. Then G has the form: 

G{p) = J2Sa{M,*p)SU 

where the sum runs over all possible energy differences a. Here denotes 
partial shifts and positive matrices. 

The map p i— > {M^- * p) preserves the energy of every state and destroys 
to some extent the coherent superpositions between them. This is shown 
by the following two extreme cases: M^- = 1 where 1 denotes the identity 
matrix. Then we have complete dephasing and obtain a mixture of energy 
eigenstates. If has only 1 as entries it is the trivial channel which does not 
affect the state at all. The map G consists of decohering channels followed 
by different energy shifts. 

To see the relation of our decomposition to the transfer function note that 
all rank-one operators |a;-|-cr)(a;| span the eigenspace of .] with eigenvalue 
(7. A time-covariant operation has to leave this subspace invariant. There- 
fore one could decompose G into a direct sum of linear maps acting on 
these eigenspaces. However, in this decomposition one would have restric- 
tions on the family B„ in order to yield a completely positive map. Therefore 
we have preferred a decomposition where each component is completely pos- 
itive. However, the Hadamard multiplication for each frequency pair {ui^ui') 
reminds to the multiplication with the value of the transfer function for each 
frequency uj. 

To generalize the decomposition to infinite dimensions we will need several 
vector- valued measures. Therefore we recall the precise definition [12]: 
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Definition 1 (Vector- Valued Measures and Densities) 

1. Let be a set and M the a-algebra of measurable subsets ofVl. For a 
Banach space B a mapping ^ : M ^ B is called a vector measure if 

7(Ujmj) = ^l{mj) 
j 

for all finite collections of mutually disjoint sets rrij . It is called count- 
ably additive when the same holds for countable sums. 

2. Let u be a m,easure on the measure space (fi, M) . A measurable func- 
tion f : Q ^ B is the Radon-Nikodym derivative of 7 with respect to v 
if 

fdu = 7(m) Vm G M . 
We also say that f is the (vector-valued) density of^ with respect to v. 

For trace-class operator- valued measures we will have countable additivity 
in the weak sense, i.e. that the scalar measure that is given by evaluation on 
observables is countably additive. For CP-map valued measures we demand 
countable additivity only after applying the maps to states and evaluating 
them on observables. 

In order to construct the infinite dimensional analogue of G^r we define a 
function 

fKAt) tr{KG{pe-'"')e'"') . 

for every observable K and state p. First consider the simple case that 
G — ApA^ where A is an operator satisfying at{A) — e~'^'^A. Then we have 

fKAt) = e'^'fKA^) . 

If G is defined by several Kraus operators Aj which are eigenvectors with 
different eigenvalues the function f^p would consist of harmonic functions 
with all these frequencies Uj. Even though we do not expect in the infinite 
dimensional case that we have Kraus operators which are eigenvectors of at it 
will turn out that we can nevertheless construct a decomposition of G based 
on the Fourier transform of fK,p- We show that for each positive bounded 
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operator K and density operator p the function fK,p is positive semidefinite, 
i.e., it satisfies 

^XkXif{tk -ti)>0 

k,l 

for all vectors x e C"* and m-tuples ti, . . . ,tm & M. with arbitrary m: 

j 

j 

where we have used 

j j 

which is just another version of the time-invariance condition. We have 
therefore 

J2^kXif{tk-ti) = 5;]ir(X(J]a;,e^^*M,e-^*'=)p(5]xK3'^M]e-^*0) 

k,l j k I 

j 

with the abbreviation 

k 

Since fK,p is positive semidefinite it defines a unique positive scalar measure 
vk,p on R by Bochner's theorem [13]. For every measurable set m C M the 
map K I— > i'K,p{m) is a positive linear functional. It is norm-continuous 
(with respect to the operator norm) since the norm of a positive functional 
is given by its value on the identity [14], i.e., by z/i,p(m). Now we restrict 
the functional to the set of compact operators where every norm-continuous 
functional is given by a trace-class operator [15]. Therefore we may define a 
positive trace-class operator p^ by 

tr{pmK) = i^K,p{m) . 
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Hence the map Gm defined by Gm{p) := Pm is a positive map on the trace- 
class operators. One can easily check that Cm is also completely positive: 
one substitutes G and e~*^* by an arbitrary tensor product extension G®id 
by and e~*^* ® 1 and considers p and K as operators on the extended space. 
Then it is obvious that the positivity argument above works similarly. We 
have Gr = G due to 

tr{Kp^) = VK,p{^) = = tr{Kp) . 

Since the trace of each G{p) is at most tr{G{p)) = tr{p) = 1, each map G^ is 
a bounded operator on T, the set of trace-class operators. Since Gk is trace- 
preserving the map m i— > Gm is formally an instrument in the sense of Davies 
(see [16], chapter 4). Think of Gm{p)/tr{Gm{,p)) as the post-measurement 
state given the knowledge that the measurement outcome a is in m. Then 
p I— > G{p) is formally the efi:ect of the measurement if the measured outcome 
is completely ignored. Even though this interpretation refers to a virtual 
measurement the virtual result a has an observable interpretation: Assume 
one observes the energy of the environment before and after it has inter- 
acted with the system. Due to energy conservation a is the energy loss of 
the environment. Note that the initial energy of the environment can be 
observed without disturbing its state since we have assumed that it is in a 
stationary state. ^ We shall call the probability measure Vp^x the "energy shift 
probabilities" in the state p. 

Each map Gm is also time-covariant: The obvious equation 

fK,p fat{K),at(p) 

implies that it is irrelevant for the measure z/^- p if at is applied to p and a^t 
to the output state (which is equivalent to applying at to K). 
We summarize the results: 

Theorem 2 For each time-covariant channel G there is a unique CP-map- 
valued-measure 

Gm 

such that 

^The idea to observe the environment in order to have a less mixed output state (which 
increases information capacities) can already be found in [17, 18]. 
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1. Each Gm is a time-covariant bounded operator on T. 

2. The Fourier transform of m ^ Gm is the function t ^ Gt with 

in the sense that 

t ^ triKGtip)) 

is for all K > and all states p the Fourier transform of the (non- 
negative) measure 

m ^ tr{KGm{p)) ■ 

The positive-map- valued measure defines a positive operator-valued-measure 
(POVM) [16] m^Qra^y 

tr{QmP) tr{Gm{p)) ■ 

In contrast to CP-map valued measures, POVMs describe only the probabili- 
ties for the measurement outcome without referring to the post-measurement 
state. The Fourier transform of this POVM will play a crucial role later. 

In order to have a stronger analogy to Theorem 1 we would like to write 
G as an integral 

G = I G,dv{a) 

with an appropriate measure v. Then the function a h— > G^ would be the 
CP-map valued Radon-Nikodym derivative of the measure m Gm, i.e., 
G(j would be the density of the measure with respect to the measure v. 
The general problem of Radon-Nikodym derivatives of instruments (CP- map 
valued measures) has already been considered in [19]. It is shown that it 
exists in the following sense: 

There is a cr-finite measure i/, a dense domain V & 7i, and a countable 
family of functions 

defined for i/-almost all a such that are linear operators V ^ H (not 
necessarily closable) such that 

k 
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and 



tr 



{KGM)= [ Y.{Al^\KAl^)duia) ^ e V 



Note that could formally be considered as the Kraus operators of a CP- 
map Ga with the "only" difference that Kraus operators are not only closable 
but even bounded. 

Equation (7) states implicitly that 



for i/-almost all a. Since it even converges for K — 1 expression (8) defines 
a bounded positive functional on the operators K. We can find a unique 
positive trace-class operator p^. such that 



for all compact K. The proof in [19] states furthermore that V can be the 
finite span of any orthogonal system of 7i. Summarizing these results, we 
have: 

Theorem 3 (Radon-Nikodym derivative of the instrument) 

For every orthonormal system {\xj)) ofH there is a family of CP-maps 
defined on the finite span of all rank-one operators \xi){xj\ such that 



for all p in the domain of Gcr and an appropriate a-finite measure v. 
However, the domain of the maps G^ can be extended: 

Lemma 1 (Extended Domains) 

Let p eT he arbitrary and v ne an arbitrary measure. If the trace-class 
operators G(j{p) are consistently defined (in the sense that they define the 
density of the measure m h- >• Gm{p) ) then we have: 




(8) 



k 
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1. The domain of can consistently be extended to pe 



-iHt 



by 



Consistency means that it is the density of the measure m i— Gm{pe ^^^) 
Similarly, we may define G„{exp{iHt)p) :— exp{—iat) exp{iHt)Ga-{p). 

2. Let I he a measurable subset ofM. and Pi be the projection onto the space 
of functions vanishing on the complement of I. Then one may extent 
the domain of G^ consistently to pPi by setting 



G^ipPi) := G,{p)Pi+, 



a -I 



and similarly GciPip) '■— Pi+aGa{p)- 



Proof: The Fourier transform of the measure 



m ^ tr{KGm{pe-'^^)) 



is given by 



tr{KGr{pe-'''')) = tr{KGt+r{p) e"^^*) 



(9) 
(10) 




The last equahty holds since Gr+t is the Fourier transform of m h- > 
evaluated at r + 1 Set p^- '■— Gaip)- We write (10) as 



m 




which is the Fourier transform of a measure with density 



a ^ tr{Kp^ e-'"*^^) . 



This proves statement (1). 

To prove (2) we have to show that 




(11) 
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holds for all K and measurable m. Consider the scalar complex- valued mea- 
sure 

I ^ tr{KG^{pPi)) . 

Its Fourier transform is 

f{t) ■.^tT{KGm{pe''''))= [ tr{Kp,e-''''e''ndi^{a), 

J m 

where the last equahty is due to statement (1). By 

the last expression is the Fourier transform of the measure 

1 G,{p)Pl+adv{<T). 

J m 

This proves that both sides of eq. (11) coincide. □ 

We would hke to characterize the maps G^r more explicitly. We already 
have done this in Theorem 1 for the finite dimensional case when the spec- 
trum of H is non-degenerate. Now we assume that H. is the set of square 
integrable functions on M 

rt := L2(M,;u), 

where p is an arbitrary measure on M defined on the Lebesgue measurable 
sets. To formalize the assumption of non-degenerate spectrum we define H 
as the multiplication operator 

Note that the Hadamard product can be generalized to infinite dimensions 
since every density operator p has a representation p — J2jPj\^j){^j\ with 
eigenvectors Then we may interpret 

P{lJ,Lj') := ^Pj(j)j{Lj)'^j{Lj') 
3 

as the entries of p. This is also the representation of p by its integral kernel, 
i.e., 

= j p{uj,uj')ijj{uj')dp{u;') . 

Accordingly, we define: 
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Definition 2 (Hadamard product in infinite dimensions) 

For two trace- class operators p, p' and a function M : M x M — >• C we write 

p' ^ M * p 

p'{uj,uj') = M{uj,uj')p{uj,uj') . 

To generalize the concept of partial shifts to L^(R, /i) is less straight- 
forward. Since they should be partial isometries, the problem is two-fold: 
Assume first that p is given by a non-constant and non-zero density with 
respect to the Lebesgue measure. Then it is intuitively clear that the gen- 
eralized "shift" rescales the function to compensate the different densities 
at different points in order to be isometric. Assume secondly that H has 
discrete and continuous spectrum, i.e., p consists of singular and absolutely 
continuous parts with respect to the Lebesgue measure. Then one cannot 
expect that the discrete part of the wave function can be shifted to the con- 
tinuous part and vice versa. The part of the wave function where this is the 
case has to be mapped to zero. 

Despite of these difficulties, the following lemma shows that our concept 
makes sense even in the general situation: 

Lemma 2 (Partial shifts) For every o" e M there is a unique partial isom- 
etry S'^. with the following properties: 

1. For all f we have 

{S^f){uj) := s^{(jj - a)f{uj - a) , 
where is an appropriate measurable "scaling" function. 

2. Under all partial isometries S„ which have a representation as in (1) 
with an appropriate scaling function the operator has minimal 
kernel, i.e., its kernel is contained in the kernel of all those S^. 

Proof: Define the translated measure Paijn) := p{m + a). Then there is 
a unique decomposition 

p^ Pc + Ps 
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where /ic is absolutely continuous with respect to jj,^ and /i^ is orthogonal 
to /Xct [20]. Here, orthogonality means that there is a measurable set B such 
that HsiB) — and HaiB) — where B denotes the complement of B. Let 
p be the density of fic with respect to /ig-. Set ■= \/p- First we show that 
Sa is an isometry of the subspace Tic given by all functions / with /(a;) = 
for LU & B: 

j \s,{u-a)f{u-a)\^dii{u) = j \s,{u)\^\f{u)\''dii,{u) (12) 
\f{u)M^)d^x,{u) = / \f{u)\^d^i, 



This integral is equal to J \f{u})\'^dii if / vanishes on B. On the other hand, 
it is zero if / vanishes on B. This shows that the space can be decomposed 
into the kernel of and a subspace where is an isometry. Hence we 
have ||>S'o-/|| < ||/|| on the whole space. Note that the equations showing 
this inequality show furthermore that S^r is well-defined with respect to the 
equivalence classes in L^(M, /i). 

Let Sa be another isometry satisfying condition (1) with the scaling func- 
tion So-. Set D :— s~^{0). Clearly, condition (1) imphes [H,Sa] — crSa and 
[H, §1] = —crSl. Hence the initial projection SlSa- commutes with H and is 
therefore the multiplication operator with the characteristic function xd of 
a measurable set D. Let / vanish on D. Since preserves its norm we have 

J \f{u)\'dix{u)^ j \f{u)rsl{u)dy^^{u). 

Since this holds for all functions which vanish on D yU-almost everywhere it 
shows that the restriction of fi to D is absolutely continuous with respect to 
lia (with the density s^). Hence the singular part of /i (with respect to iia) 
vanishes on D and we have 



J \f{u;)\'df,{u;)^ J \f{u)fdi,,{u;). 



As can be seen in eq. (12), this shows that the norm of / is also preserved 
by Sa- Hence the isometric subspace of Sa is a subspace of the isometric 
subspace of Sa and the kernel of Sa is a subspace of the kernel of S'o-.D 

Fortunately, these shifts have the following property: 
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Lemma 3 The partial shift of Lemma 2 satisfy 

Sl — S^a- ■ 

Proof: Rewriting the inner product of Ti as an integral one checks easily 
that Sl has to be a translation by —a with an appropriate scaling function. It 
is easy to see that S-^Sa = SlS^: On the kernel of S'o- both products coincide 
trivially. On its orthogonal complement, i.e., the isometric subscape of the 
operator SlSa- is the identity. But this must also be true for S-^S^ because 
the isometric subspace of S-a must contain the isometric subspace of since 
S-a- has minimal kernel in the set of all partially isometric translations (see 
property (2) in Lemma 2). Using general properties of partial isometries the 
image of satisfies therefore 

ImSl = Im(SlSa) = Im(S^^Sa) < ImS^^ . 
Reversing this inequality by taking the orthogonal complement we obtain 

kerS(j > kerS'L^ . 
Because this holds for all cr e M we have 

kerS-cr > kerSl . 
Since S^^ has minimal kernel we conclude 

kerS-cr — kerSl . 

Due to the uniqueness of the maximal element in the set of all partially 
isometric translations we have 

□ 

Note that we can only write G^ip) in the form S^iM^ * p)Sl if the density 
operator G„{p) acts only on the image of So- (which coincides with the image 
of the projection S^Sl). But this is indeed the case: 

Lemma 4 If p is in the domain of its output p^ '■— Ga{p) satisfies 

SaSl-PaSaSl- = Pa ■ 
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Proof: We will show Ga{p)SaSl = Ga{p)- The statement SaSlG^ip) = 
Gcr{p) follows similarly. Due to Lemma 3 we have S^^Sl = S^_^S_(^, which is 
the projection onto the isometric subspace of S^^-- It can explictly be given 
as follows (see proof of Lemma 2): Let pi := p^c® l^s be the decomposition 
into the absolutely continuous and singular part of pi with respect to /i-o-- 
Let C be a set such that pi_^{G) = and pis{G) = 0. Then S'^_„S-cr = Pc- 
Due to /i(C — a) = p,-a{G) = we have Pc_^ = 0. This impUcs 

G,{p)Pc^G,{pPc_,) = Q. 

□ 



For the construction in the sequel we choose the ONS in Theorem 3 such 
that its finite span contains a vector \ip) with ■ip{u!) ^ for /x-almost all a;. 
Now we can define a function 

: M X M ^ C 

by 

I A k(uj,u') 

where the function k represents the trace-class operator Slp^Sa- Even though 
this function is not well-defined, the calculations below show that the freedom 
of choosing is irrelevant. 

By construction and using Lemma 4, we have 

G^{p)=p^ = S„{M,*p,)Sl. 

We conclude 

G(\^)m)^ J S,(M,*(\^l;){^l;\))Sidu(a) 
with l'^) as above. Using Lemma 1, part (1), we have 

for all s, t G M. The finite span of these rank-one operators is dense in T since 
H is the multiplication with the identity and the span of all e'xp{—iHt)\ilj) 
is therefore dense in 7i. We conclude: 
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Theorem 4 (Dephasing - Energy Shift Representation) 

There is a family of functions {M„)aeM 

: M X M ^ C , 

and a a-finite measure v on'M. such that for a dense set of density operators 
p the decompostion 

G{p) = J S„{M,*p)Sldu{a) 

holds. 



We would like to have an analogy to the statement that M„ is in the finite 
dimensional case a positive matrix. However, we have specified only as a 
function and not as an operator. It seems straightforward to consider it as 
an operator by 

(M,V')H := j M,{uj,J)'il^{J)dpi{J). 

But it is easy to see that this is not in general well-defined. Let, for instance, 
G be the identity on the density operators on /^(Z). Then the decomposition 
of G reduces to G{p) — Mq * p with Mo(uj,u;') — 1. The formal matrix 
multiplication of this "all-one" matrix with any square-integrable sequence 
leads to infinite coefficients. However, formally it is like a positive operator 
in the following sense. 

is by construction the trace-class operator p„. Translated to finite dimensions 
this expression to 

where D is a non-singular diagonal matrix (in analogy to the statement that 
vanishes almost nowhere. Such a matrix can only be positive if Mg. is 
itself positive. In this sense we consider as positive even though it is not 
an operator. In analogy to the finite dimensional case, an exact measurement 
of the energy of the environment (before and after the interaction has taken 
place) would lead to a map 

S,{M,*p)Sl, 
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which is only a dcphasing channel up to the known energy shift. 

It should be noted that (depending on u) the integral may represent a 
finite, or countable infinite sum, or a continuous integral even though the 
Kraus representation of G is always possible with a countable sum [2]. 

7 Superpositions between states in the same 
orbit 

The "Hadamard-channel" p M*p with a positive matrix M allows perfect 
classical information transfer by using energy eigenstates as logical states. Its 
quantum capacity depends on the destruction of off-diagonal elements, i.e., 
whether the output states are more or less stationary states. Explicitly one 
has the following lower bound on the quantum capacity: 

Lemma 5 In dimension n, the quantum capacity Q of the channel p i-^ M*p 
with M{u,u) = 1 for all uj E fl satisfies 

Q>\og^{n)-S{M/n) 

where S{.) denotes the von-Neumann entropy. 

Proof: Let 

(13) 

i 

be a decomposition into mutually orthogonal rank-one operators with non- 
normalized vectors \mj). Let Dj be the diagonal matrices having the coeffi- 
cients of \mj) as entries. Then we have 

G{p) = ^oPD] . 

3 

The quantum capacity is given by the maximum of the coherent information 
[23]. Let 10) (01 a state on C"®C" and p its restriction to the right component. 
Then the coherent information is defined as 

S{G{p))-S{id®G){\(t>)m, 
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where S denotes the von-Neumann entropy. Consider the maximally entan- 
gled state 

Note that p is the maximally mixed state and that it is preserved by the 
channel. Hence we have G{S{p)) — log2(n) and 

{id ® G)(|0)(0|) = ® D,)\<p){<p\{l ® D}) . (14) 
j 

Furthermore we have 

(0|(1®dJ)(1® =tr(DjA) = {m,\mi) = 0. 

Therefore the sum (14) is already the spectral decomposition of {id®G){\(f)) (01) 
The eigenvalues are the square of the length of each vector (1 Dj)\(f)). It 
is given by {mj\mj) /n. Hence the entropy of the output state is the von- 
Neumann entropy S{M/n). □ 

Here the destruction of quantum information coincides with the destruc- 
tion of timing information in the sense that states on the same orbit of 
become less distinguishable. It is obvious that decoherence with respect to 
the energy eigenbasis leads destroys always timing information in this sense. 
Here we want to address the question what happens when one considers or- 
thogonal states on the same orbit of the time evolution as reference basis. 
Is it possible that a timc-covariant channel destroys superpositions between 
them without affecting the basis states? Consider two vector states |0o) and 
\4>s) := exp(-iifs)|0o)- Clearly, G'(|0o)(0o|) = |0o)(0o| implies 

^.(M.*(|0o)(0o|))5^ = O 

for all o" 7^ 0. This follows from the observation that the output state can 
only be pure if all output operators are linearly dependent. Every cr 7^ 
would lead to a different output state. Furthermore we have 

Mo * (|0o)(0.|) = Mo * {\ct>o){My' = |0o)(0.| ■ 

One can conclude easily that every density operator acting on the subspace 
spanned by |0o) and l^^) is preserved by G. A little bit more general, one 
has: 
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Theorem 5 Let G be a time-covariant CP-map on a system with 
non- degenerated Hamiltonian. If there is any pure state \4>){(p\ with 

then G leaves all density operators invariant which act on the Hilbert suhspace 
spanned by all exp{—iHt)\(f)) . 

The essential argument above is that there exist no states which are 
invariant with respect to S„ for some a ^ 0. An analogue statement would 
not be true for covariance with respect to a discretized time evolution: 

Consider the Hilbert space of square summable two-sided sequences 
and the discrete time evolution given by the translation {Uip){n) -0(71 + 1). 
Let (|e„))„gz be the canaonical basis of /^(Z). Then the channel 

|e„)(e„|p|e„)(e„| 

leaves all basis states \en) invariant but destroys all superpositions. 

The assumption that a pure state is preserved by G is rather strong. 
Actually we want to figure out whether distiguishability of different states 
in the same orbit may be conserved even though their superpositions are 
distroyed. We consider the following extreme case: 

Definition 3 (Reliable timing) 

A time-covariant channel G has the "reliable timing property" (with respect 
to the time s) if there exists an input density operator p and a real num- 
ber s such that for ps := Ois{p) the outputs G{p) and G{ps) are perfectly 
distinguishable, i.e., the density matrices are mutually orthogonal. 

It is easy to verify that the input state p can be chosen to be pure. To 
justify the definition, we show that this property appears in the following 
situation: Assume a sender, say Alice, wants to send a signal to a receiver, 
say Bob. Assume furthermore that it should be guaranteed that Bob receives 
the signal in the time interval [ti,t2]- This requires that the physical state 
p of the signal (when it is sent) is perfectly distinguishable from the time 
evolved state ps with s := ti — t2. This is due to the fact that the following 
"measurement" distinguishes between them: wait the time ti and ask Bob 
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whether he has ah^cady received the signal. If the medium between Ahce and 
Bob modifies the signal, we may model this by a time-covariant operation G 
which has clearly to preserve the distinguishability of the states p and Ps- 

It is clear that reliable transfer of classical information requires two states 
p and 7 such that G{p) and (^(7) are mutually orthogonal. The remark above 
shows that there are situations in classical information processing where two 
orthogonal output states should exists on the same orbit since reliable timing 
requires this feature. To see that non-trivial channels with this property exist 
one may construct a channel of the form 

j 

where the aj are chosen such that {(f)\S'l.Srj <i>) = for i 7^ j. This condition 
is, for instance, satisfied if the minimal distance between the values aj exceeds 
the spectral bandwidth of the input state. Real physical channels will satisfy 
condition 3 at most approximatively. 

In order to consider the quantum capacity of reliable timing channels we 
choose p and ps = UgpUj with Ug := exp(— iiJs) such that p and ps are 
mutually orthogonal. For simplicity we assume that the dynamical evolution 
of G{p) is periodic. We obtain a set 

G{po), G{ps), G{p2s), G{ps(N-l)) 

of N mutually perfectly distinguishable density matrices. Let P be the sup- 
port of i.e. the smallest projection with PG{p) — G{p). Then the 
projections UgjPUlj with j = 0, . . . , — 1 are mutually orthogonal. 
Now we consider the channel with input 

= span{\(f)sj)}j=o,...N-i ■ 

It is clear that there exist pure input states \ipo), . . . , lipsiN-i)) for which the 
orthogonality of the output states is satisfied. By sending either of these 
N states G can transfer logsA^ bits of classical information. We want to 
derive sufficient conditions under which G allows also to send superpositions 
between the chosen basis states. For doing so, we restrict G to a channel K 
on N X A'"-density matrices as follows. 
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The input is restricted to the span of all \4>sj) with j = 0, . . . , A?" — 1. The 
corresponding output space is 

We may consider this space as the tensor product 

nr = C^ PH 

if we identify the spaces UsjPUlpi, for j 7^ with PH. via arbitrary unitaries. 

It is straightforward and convenient to chose the isomorphisms Usj with j = 
0, . . . ,N — 1. For any output density operator 7 acting on Tir the entry 
corresponding to \j){k\ of its restriction to is given by 

tr{U,kPUlj^) . 

To determine the channel K we have to compute all values 

triU,kPUl^G{\<Psi){<f>sm\)), 

for k,j, 1,171 e {0, . . . , N — 1}. Due to the reliable timing property we know 
that each state \(f>sj) {(f>sj\ leads to the output state Furthermore the 

operator \4>si){4>sm\ leads with certainty to a multiple of Roughly 
speaking, the reason is that a CP-map which maps the states \m){m\ and 
|/)(/| onto itself maps also onto multiples of itself. This could, for 

instance, be shown by reformulating the CP-map as the restriction of an 
appropriate unitary map (see [14]). Hence we have only to determine which 
factors the off-diagonal terms obtain. Due to the symmetry of the channel 
with respect to time translations Usj we have only to evaluate 

v{j) := tr{UsjPG{\<f>o){<f>sj\)) = triPG{\<l>o){<f>sj\Usj)) ■ 
Then K is explicitly given by 

K(^Cjk\j){k\) ^^Cjkv{j - k)\j){k\ , 

jk jk 

This shows that the channel is a "Hadamard channel" even though we have 
not chosen the energy states as reference basis but the states {(psj) instead. 
Furthermore, K is given by Hadamard multiplication with a circulant matrix 
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V with entries Vij := v{i — j). In order to calculate v we use the explicit 
form of G according to Theorem 3. Using Lemma 1 this yields 

v{j) = I tr{GMo){^o\)e-"''' du{a) = f^Mi^A-^j) , 

i.e., the Fourier transform of the energy shift probability measure z^|<^o)((^(,|,i 
evaluated at the points —sj. The eigenvalues qo, . . . , Qm-i of V/n are given 
by the inverse Fourier transform of v: 

j=0 

Using Lemma 5 we conclude: 

Theorem 6 (Q-capacity of channels with reliable timing) 

Let G be a time-covariant channel with reliable timing property with respect 
to the time s and the initial state |0o)- Let the dynamics be periodic with 
period length sN . Then the quantum capacity can only be zero if the Fourier 
transform p of the probability measure p{m) := tr(G'„t(|0o)(0o|)) has zeros at 
sj for j = 1, . . . , N — 1. Otherwise the quantum capacity is at least 

log2-/V- log2 qk,, 

k 

where q — (go, ■ • • , qN-i) is the discrete Fourier transform of the evaluation 
of f\4,o){MA ^"^ ^ points 0, -s, -2s, ...,-{N- l)s. 

Of course there are many possibilities to define finite dimensional channels 
from the original one. However, to study whether G destroys superpositions 
between different states in the orbit ("different pointer states of a clock") 
one has always to consider 

G(|0o)(0.,-|) = G-.,-(|0o)(0o|)e^''^^ 

This shows that the Fourier transform of the measure {Gm) is decisive. For 
deriving the lower bound above we have evaluated it for the observable 1 
which leads to the Fourier transform of the POVM {Qm) defined in the last 
section. However, 1 is not necessarily optimal for detecting superpositions 
between the output states. 
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8 Gaussian channel on a single mode 



The Hilbert space of a single mode in quantum optics is i^(No), the set of 

square summable sequences. It can also be interpreted as the energy levels of 
a harmonic oscillator. A time-covariant channel which is often considered is 
the following. Using position and momentum observables X and P, respec- 
tively, we define a translation (01,02) G IR^ by the unitary transformation 
exp(i(aiX + a2-P)), where X :— a''' + a and P :— {a'' — a)/i. It is conve- 
nient to introduce variables r G and z on the complex unit circle F by 
ai + ia2 = rz. Rewriting the translation with operators and a and the 
parameters z,r we obtain the translation operator (compare [24]) by 

D{z,r) := exp(r(^a^ - za)) = g-VSgr-Wg-r^a _ 

Since the global phase factor is irrelevant we will, in abuse of notation, use 
D{z,r) for the term without this factor. 

Let the gaussian channel be given by a random displacement rz accord- 
ing to the two-dimensional rotation invariant gauss distribution. Note that 
rotation symmetry is necessary in order to obtain a time-covariant channel 
since the Hamiltonian H := diag{0, 1,2... .) corresponds to a rotation in the 
"phase space" . Note that the possible energy shifts are in Z, hence we expect 
a countable sum of which are (in contrary to the general case) not only 
densely defined. 

The whole channel is given by (see [25]) 



D{z, r)pD{z, r) dzp{r)dr 



with p{r) := exp(— r^/(2s^))r/s^ and s denotes the standard deviation. 

As usual, we introduce the creation operator a'' and the annihilation 
operator a by a\j) := V7b ~ 1) ^or j > I and a|0) = 0. Then we write 
D{z,r) as the power series 

m>0 ■ n>0 

It decomposes canonically into terms Dfj{z,r) with cr G Z satisfying the 
commutation relation [H, Du{z,r)] = aD„{z,r) if we define 



V (rzJ-"(at)-+-(-r^)'-a- _ ^ r-l^i^^^fr^r 



n>— <T,0 ^ ' n>— cr.O 
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Note that Dcr{z, r) maps states with energy j onto states with energy j + a. 
This suggests already that they may correspond to the maps of Section 
6. The operator 

is independent of z. We conclude therefore that terms of the form D^{z, r)pD^/{z, r) 
cancel for a ^ a' after integration over all 2; e F. We obtain 



Consider the case cr > first. Then we have 

1 (j + a)! 



for all n < j. For the case cr < we obtain 



For (J > we conclude 



where 



' ^ (j-n)!(n + (7)!n!' 

is a Laguerre polynomial [26]. 
Hence we have 

J + 1) . . . + (7) 0' + 1) . . . if + a) ' 

For negative a we may express similarly by Laguerre polynomials together 
with a different factor. Hence the decomposition of the single mode gaussian 
channel can be given in a closed form even though it leads to a less familiar 
representation. 
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9 Conclusions 



We have shown that every time-covariant CP-map has a representation as 
an integral over a family of CP-maps. If the spectrum of the system Hamil- 
tonian is non-degenerated, each of these components consists of a Hadamard 
multiplication with a positive operator followed by an energy shift. Formally, 
the output of the channel is an unselected post-measurement state. The mea- 
sured (and ignored) quantity is the energy which has been transferred to or 
from the environment. Conditional to the measurement outcome a different 
dephasing channel is applied. 

Furthermore we have addressed the question to what extent covariant 
channels can destroy superpositions between N mutually orthogonal states 
on the same orbit of the time evolution. For a specific type of channels (with 
"reliable timing property") we have shown that the general decomposition 
helps to derive lower bounds on the quantum capacity. 

The decomposition presented here may be a helpful approach to describe 
a rather general type of decoherence and relaxation phenomena. 

We have calculated the decomposition explicitly for a rotation invariant 
gaussian channel acting on the state space of a single mode Fock space. The 
dephasing operations can then be described using Laguerre polynomials. 
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